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Abstract: The hadronic decay rate of the r lepton serves as one of the most precise 
determinations of the QCD coupling a s . The dominant theoretical source of uncertainty at 
present resides in the seeming disparity of two approaches to improving the perturbative 
expansion with the help of the renormalisation group, namely fixed-order and contour- 
improved perturbation theory. In this work it is demonstrated that in fact both approaches 
yield compatible results. However, the fixed-order series is found to oscillate around the 
contour-improved result with an oscillation frequency of approximately six perturbative 
orders, approaching it until about the 30th order, after which the expansion reveals its 
asymptotic nature. Additionally, the renormalisation scale and scheme dependencies of 
the perturbative series for the r hadronic width are investigated in detail. 
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1. Introduction 

Already more than a decade ago it was realised that the hadronic decay of the r lepton 
could serve as an ideal system to study low-energy QCD under rather clean conditions [1-4]. 
In the following years, detailed investigations of the r hadronic width as well as invariant 
mass distributions have served to determine the QCD coupling a s to a precision competitive 
with the current world average [5-10]. The experimental separation of the Cabibbo-allowed 
decays and Cabibbo-suppressed modes into strange particles [5,11,12] opened a means to 
also determine the quark-mixing matrix element \V US \ [13-15] as well as the mass of the 
strange quark [16-23], additional fundamental parameters within the Standard Model. 

The starting point for a QCD analysis of the t hadronic width R T is the finite energy 
sum rule (FESR) [1,24] 

so 

T(r — > hadrons v T ) _ n f . . . . —1 



= E^(sq) = j w T (s) p T (s) ds = — ^ j> w T (s)U T (s)ds . (1.1) 



T(r ->• HV^Vt) 

| s |= ao 

The FESR can easily be derived from Cauchy's theorem and the fact that the exact, non- 
perturbative correlation function IT r (s) is analytic in the whole complex s-plane, except for 
the positive real axis, where it may have poles and cuts. The correlator II T (s) of hadronic 
QCD currents receives contributions from vector and axialvector correlation functions for 
both {ud) as well as {us) flavour content, dressed with the appropriate quark mixing matrix 
factors: 



n T ( s ) = |K<^ TiXM + K d (s) +\Vus\ z Ks(s) + Ks(s) • (1.2) 



rV , n \ , rrA 
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Furthermore, p T (s) = ImIl T (s + iQ)/ir is the so-called spectral function which contains all 
the physical information. 

In the case of hadronic r decays, s is the invariant mass of the final hadron system, 
w T (s) is the relevant kinematic weight function, 

„*(.)_ - -iV(l + 2 ±), (1 . 3) 

so V s oJ V s oJ 

and so = ^ T 2 i though for the moment the somewhat more general expression will be 
kept. The fortuitous fact about R T is that the weight function introduces a double zero at 
s = so, such that there is no need to evaluate II T (s) close to the real axis where perturbation 
theory becomes problematic. In writing eq. (b_l), a small contribution from longitudinal 
correlators has been omitted. Being suppressed by the light quark masses, it is of no 
relevance for the purpose of this work which will concentrate on the massless correlators, 
but it will be reconsidered in a subsequent publication. 

Since the hadronic r decay rate is a physical, measurable quantity, the corresponding 
perturbative QCD expression can be improved with the help of the renormalisation group 
equation (RGE). In the course of calculating i?^(so), one also has to perform the contour 
integration over the circle with radius so in the complex s-plane. Then the question arises 
in which order both operations should be performed. First calculating the contour integral 
and then performing the renormalisation group improvement goes under the name of fixed- 
order perturbation theory (FOPT) whereas the second approach of first resumming the 
expansion and afterwards integrating over the contour has been suggested in [25,26] and is 
termed contour-improved perturbation theory (CIPT). Both approaches will be investigated 
in much detail below. 

Numerically, it is found that the CI perturbative expansion displays a better behaviour 
than the corresponding FO series [26]. In particular, employing the known perturbative 
results up to order as well as estimates for the contributions at order a* and a^, 
it is found that the difference between the two approaches is much larger than would 
be naively expected, based on the last included terms in the expansion. In the past, 
this apparent ambiguity has cast doubts as to which way of resumming the perturbative 
series is preferable and it represents the dominant source of theoretical uncertainty in the 
determination of a s from the hadronic t decay rate up to this day. 

Below, it is demonstrated that in fact both approaches to the renormalisation group 
improvement of the perturbative series yield compatible results. However, the fixed-order 
series is found to oscillate around the contour-improved result with an oscillation frequency 
of approximately six perturbative orders, approaching it until about the 30th order, after 
which the expansion reveals its asymptotic nature. To this end, in the next section, the 
perturbative expansion of the vector correlation function will be reviewed, and in section 3 
the two approaches of improving the perturbative series in the r-decay finite energy sum 
rule will be discussed. In section 4, a simplified example will be analysed, in which only the 
first coefficient of the /3-function /?i is non-vanishing. Whereas this case can be treated fully 
analytically, it shares all the essential features of the complete QCD expression. The more 
realistic case which includes all currently known terms in the /3-function is then treated in 
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section 5. In section 6 and 7, the influence of variations of the renormalisation scale as well 
as the renormalisation scheme on the hadronic r decay rate are investigated, and, finally, 
in the conclusions my main findings of this work will be summarised. 



2. Perturbative correlator 

By far the most important contribution to the hadronic r decay rate is due to the purely 
perturbative vector (axialvector) correlation function in the limit of vanishing quark masses 
[1]. In this limit vector and axialvector correlators yield identical contributions. Therefore, 
in what follows, only the vector contribution will be investigated. A discussion of quark 
mass corrections shall be presented in a forthcoming publication, but for the remainder of 
this work, the quark masses are assumed to be zero. 

In the massless case, the perturbative expansion of the vector correlation function II(s) 
exhibits the general structure: 1 

pj oo n+1 

n(a) = - ^ E a l E c ^ L " where L = ln ^ > 

77 n=0 k=0 ^ 

and a„ = a s (p 2 )/ir with p being the renormalisation scale. (The global minus sign has 
been introduced for later convenience.) II(s) itself is not a physical quantity. However, this 
is the case for the spectral function p(s) as well as for the Adler function D(s): 

, oo n+1 

D(s) = - s ± U(s) = Y^J2<J2 kCnkLk ~ 1 - ( 2 - 2 ) 

n=0 k=l 

It should be clear that the coefficients c h q do not contribute to both the spectral function 
(because they are real) and to the Adler function (because of the factor k). Thus the 
coefficients c n o, which actually are found to depend on the renormalisation prescription, 
can be considered as "unphysical" in that they do not appear in measurable quantities. 
As physical quantities, both D(s) as well as p(s) have to satisfy a homogeneous RGE: 



d \d(s) 



< 




0, (2.3) 



dp [p(s)j 

where (3(a) is the QCD /3- function, defined as: 

- p ^ = (3{a) =Aa 2 + /? 2 a 3 + fta 4 + fta 5 + .... (2.4) 

Numerically, for N c = 3 and in the MS-scheme [27] the first four coefficients are given 
by [28-30]: 

a — 11 1 v a _ 51 19 at a _ 2857 5033 AT , 325 AT 2 

Pi- -2-3*1/, fh-—-j2^f' P3 ~ ST ~ W N f + T728 N f ' 

a _ 149753 , 891 /■ ( 1078361 , 1627 r \ AT _L ( 50065 , 809 /- \ AT 2 ^. 1093 AT 3 (O K\ 
~ ~yW- + 12" ^3 - 20736 + + ^20736 + T296 ^ ) ^ f + 93312 iV / ■ \ 2 - b ) 



1 Since in the rest of this work, we will only be concerned with the vector correlator, the superscript V 
will be dropped, and also the flavour content does not matter in the limit of vanishing quark masses. 
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The RGE puts constraints on the coefficients c n k- Considering the coefficients c n \ to 
be independent, all other coefficients c n k with k = 2, . . . , n + 1 can be expressed in terms 
of the c n i and 0- function coefficients. Up to order , the RG constraints lead to: 

C22 = - y en , c 33 = y| en , c 32 = - | (/? 2 en + 2/?i c 2 i) , (2.6) 

/?i 1 

C44 = ~ 32 Cl1 ' C43 = 24 ^ 2 Cn + ° 21 ^ ' ° 42 = ~ 4 ^ 3 ° n + 2 ^ 2 ° 21 + C31 ^ ' 

Furthermore, the coefficients c nj „+i = for n > 1. The independent coefficients c n i are 
known analytically up to order [31, 32] and at N c = 3 in the MS-scheme take the 
following values: 

coi = cn = 1 , C2i = ^ - HCs - (n " fCa)^/ , 

_ _ 87029 1103 /- , 275/- / 7847 262 a , 25 a "\ at , f 151 19 a \ jy-2 / 9 7 \ 

c 3i - ^88 — Cs + — Cs - - — Cs + tCsJ% + (^162 - w&) N f • (2.7) 

For the next five- and six-loop coefficients C41 and C51, estimates employing principles of 
"minimal sensitivity" (PMS) or "fastest apparent convergence" (FAC) [33,34], together 
with recently computed terms of order a^Nj, exist which read [35,36]: 

c 4 i = 27 ± 16 , c 5 i = 145 ± 100 . (2.8) 

For illustrative purposes, also the central values of these estimates will be taken into account 
in the analysis below. 

Since the Adler function D(s) satisfies a homogeneous RGE, the logarithms in eq. ( [2.2; ) 
can be resummed with the choice fj? = —s = Q 2 , leading to the simple expression: 

AT 00 



12vr 2 

ra=0 

where oq = a s (Q 2 ) /ir. As is again apparent from this equation, the only physically relevant 
coefficients are the c n \. 

3. Hadronic r decay rate 

In principle, now one could proceed by inserting the general expression ([2.1|) for I4(s) into 



the contour integral of eq. (1.1). However, it is advantageous to rewrite the FESR for the 



r decay rate in terms of the Adler function D(s) by partial integration: 2 

R?(s ) = - 6m <j> ^ (1 - xf (1 + x) D(s x) , (3.1) 



x =1 



with the new dimensionless integration variable being x = s/sq. For FOPT, order by order 
R^{s ) and R?(s ) can be shown to be identical as they should, but applying CIPT, the 



In writing eq. (|3.l[), \V u d\ + \ V US \ = 1 has been assumed, which numerically is satisfied rather well. 
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perturbative expansion for R^(sq) is less well behaved than the one for (sq). Therefore, 
only the comparison of R®(sq) in FOPT and CIPT will be investigated in this work. 

Inserting the expression ( |2.2| ) for D(s) into the contour integral of eq. ( |3.1[) , one finds: 



= ^ E «m E * ^ f d i (i - -) 3 (i + *) ^ (^) 

n=0 fe=l [x J =1 \ * ' 

E^E^^Ef , j lnfe "'" 1 5 f f a - -) 3 a + *) 



2m 

n=0 fc=l Z=0 |x|=1 

00 n+1 fc— 1 



= iv c E^E fc ^E , H^S- ( 3 - 2 ) 

n=0 fc=l Z=0 V ' ^ 

In the last line, the contour integrals J\ are defined by 



Jl = hi /^(l-^) 3 (l+^)ln'(-^) = ^ 

|x|=l 

with the required integrals I\ m being given by 



1,0 + 2/0-21,3-/; 



(.1 



(3.3) 



I l>m = i l j a 1 e tma da = i (-^) T(l + 1, -iam) 



where T(/ + 1, z) is the incomplete T-function, [n] denotes the integer part of n and m > 1. 
For the case m = 0, one obtains /^o = + (—1)'] vr' +1 /(/ + 1). Explicitly, the first few 
of the integrals J/, which are needed up to order af, read: 

7—1 7 — 19 7_ 265 1 _2 t _ 3355 , 19 Jl fo c\ 

in agreement with ref. [26]. 

As discussed above, the Adler function D(s) is a physical quantity, and thus satisfies 
a homogeneous RGE. Therefore, the logarithms in eq. ( |3.2j ) can be resummed. Let us first 
concentrate on FO perturbation theory, which amounts to taking the choice fi 2 = so in the 
last line of Q3.2| ), which then leads to 

oo n+1 

r? > FO (s ) = a+e < E k c «fc J fc-i ' ( 3 - 6 ) 

n=0 k=l 

with a SQ = a(so). Inserting the relations for the c n k together with the corresponding 
contour integrals Jk-i, up to order , this leads to: 



R? ' FO (s ) = JV C | coi + en a so + (c 2 i + ^Picu)a 2 
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+ (c 3 i + gAci + (i/3 2 + (H - fj)/3i 2 )cii)^ 
+ (c4i + f Ac 31 + (i/? 2 + (H - ^)/? 2 )c 2 i 

+ (1/33 + (ff - + (H - !-> 3 >u)<4 | (3.7) 

( 19907171 i „ 22683 Jl ( 345405 729 „2\ a , 12825 ^ \ J I /a fi \ 
+ ^ 6144 + C 41 - "W 71 ~ \-T2&- 16" J Cs + — C 5 J a S() ^ . (3.8) 

For the second step, the coefficients of the /3-function (|2.5| ) as well as the c n i-coefficients of 
eq. 

Q h ave been employed. An important remark is in order here. As should be clear 
from eq. ( |3.7|) , each order n in FOPT depends on all coefficients en with i < n. This will 
play a crucial role in what follows. Before analysing the expressions numerically, the second 
approach of resumming the logarithms, namely CI perturbation theory, will be introduced. 

For CIPT, the logarithms in eq. (|3.2| ) should be resummed before performing the 
contour integral. This can be achieved with the choice /j, 2 = —sox in the first line of (|3.2|), 
which yields 

oo 

R?> cl (s ) = N C Y^ c nl JZ(s ) . (3.9) 

n=0 

Here, the equation has been rewritten in terms of the contour integrals J^(so), defined as: 
J%(so) = ^~ I — (1 - xf (1 + x) a n (-s x) . (3.10) 

l-Kl J X 
1*1=1 

In contrast to the FO case, for CIPT each order n just depends on the corresponding 
coefficient c n \. Thus, all contributions proportional to the coefficient c n \ which in FOPT 
appear at all perturbative orders equal or greater to n are resummed into a single term. 
This is related to the fact that CIPT resums the running of the QCD coupling along 
the integration contour in the complex s-plane. To perform the contour integration, we 
have to analytically continue the strong coupling a s to the complex s-plane, but this is 
straightforward since the dependence of a s on s is only logarithmic and it has the same 
cut structure as II (s). 

4. A simple example 

Before embarking on the general integrals J°(so), as a first step, let us investigate a sim- 
plified example which can be treated fully analytically, namely the case where only the 
first coefficient of the /3-function, /?i , is non-zero, and all higher order /^-coefficients vanish 
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identically. Then the running coupling in eq. ( 3.10| ) can be expressed as 

a n (-s x) ■ 



so 



1 + \ a S0 ln(-x) 

and the integrals J^(so) can be calculated analytically with the result: 

Jn(so) = < J n fi( K ) + 2 Jn,l(K) - 2J„ i3 (k) - J n>4 (/«) 

Here, the constant k = /?ia S0 /2, and the integrals J n ,m(n) are found to be 



(4.1) 



(4.2) 



n-1 



2mKT(n) 



n-1 



Ei 



z=i 



2 = l+i7TK 



- 2 = 1— i-KK 



(4.3) 



with Ei(z) = J e z / zdz being the exponential-integral function and m > 1. Again, the case 
m = which is also required has to be treated separately and yields: 



Ji,o( K ) 



1 . (l+inn) 
in ■ 



Ittik (1 — inn) 



n = 1 , 



, 1 — i7TK 



l-n 



(1 + ITTK 



, 1—n 



n > 2. 



(4.4) 



(4.5) 



2irin(n — 1) 

The analytic expressions for the integrals J°(so) can readily be compared with their 
perturbative expansion which is obtained by expanding eq. (4.1) in powers of k = (5\a so /2, 
before performing the contour integration: 



n 00 

■*(»> = £ 

v ' 1=0 



T(n + l) 
l\ 



Jl (-«)' 



(4.6) 



Employing the asymptotic expansion for the exponential-integral function in eq. (4.3) as 
well as expanding eqs. ( ^4.3[) to (|4.5|) in powers of k, one can verify the agreement between 
both representations for J%(sq). 

Let us now turn to a numerical comparison of CI versus FO perturbation theory in 
the simplified example. To do this, an input value for a s is required and the very recent 
result a s (M^) = 0.34 [6] will be employed for a So . The CI expression of eq. Q4.2| ) and the 
perturbative expansion of the FO result of eq. (|4.6j), together with eq. (|3~9|), then lead to: 



,i 



R D f = 3 

T ! PI 



= 3 

T, Pi 



a a a a a a 

1 + 0.1455 + 0.0305 + 0.0136 + 0.0060 + 0.0031 

1 + 0.1082 + 0.0609 + 0.0254 + 0.0082 + 0.0024 



3.5961 



3.6154. 



(4.7) 



(4.8) 



From (|4.7| ) and (4.8) one observes that the CI expansion appears to display a better con- 
vergence, although the last included term of the FO expansion is somewhat smaller. Al- 
together, the difference between the two results amounts to 3 • 0.0064, about 2-3 times the 
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r 
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Q 

o 


0.108zz5 


0.060933 


A AO C OA A 

0.025394 


A AAO ICC 

U.OOoloo 


A A AO A 1 A 
0.002419 


A AAACAT 

-0.000607 


A AA A 1 AO 

-0.004193 


A AA ^ccn 

-0.004669 


9 


10 


11 


12 


13 


14 


15 


16 


-0.000104 


0.003728 


0.001815 


-0.001991 


-0.001863 


0.000850 


0.001407 


-0.000275 


17 


18 


19 


20 


21 


22 


23 


24 


-0.000939 


0.000032 


0.000588 


0.000051 


-0.000354 


-0.000066 


0.000209 


0.000057 


25 


26 


27 


28 


29 


30 


31 


32 


-0.000122 


-0.000043 


0.000070 


0.000030 


-0.000040 


-0.000020 


0.000023 


0.000013 



Table 1: Perturbative order n and corresponding terms in the expansion of R T 'a /3. 



last term in both expansions. Thus, the agreement between both approaches to performing 
the renormalisation group improvement is reasonably good. 

Nevertheless, if we were to assume that all coefficients c n \ with n > 6 vanish identically, 
the CI result of eq. fl4.7| ) would represent the complete answer. Then the question arises, 
how the FO result approaches this value at higher orders in perturbation theory. Since 
with eq. (|4.6j ) an all order expression is at our disposal, the question can be answered by 

The numerical result of this 



inserting (4J3) into eq. ( |3.9|) and reexpanding in powers of a so . 
exercise is displayed in table |l| up to the 32nd order. The entries are the perturbative order 
n at the one hand and the corresponding term in the expansion of i?^^°/3 on the other. 
Therefore, the first five entries agree with the corresponding terms in the square bracket 
of eq. (p|). 

By inspecting table [l], it is observed that after the first five positive terms, there are 
four terms with negative sign. After that, the sign changes at about every second term. 
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Figure 1: Difference (R^° - R^)/3 as a function of the perturbative order n up to which the 
terms in the FOPT series are summed in the example with only j3\ non-vanishing. 
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Generally, the size of the terms decreases with increasing order, but not from one term 
to the next. A graphical account of these findings is presented in figure [l], which shows 
the difference (R^'J® — R^^)/^ as a function of the perturbative order n up to which 
the terms in the FOPT series are summed. As can be seen from this figure, the sign 
changes lead to FO results which oscillate around the CI value. Due to the sign change 
at about every second term, the oscillation frequency is approximately four perturbative 
orders. Furthermore, the magnitude of the oscillation is tending towards zero, such that 
the FO value approaches the CI result at large orders. Incidentally, the maximum of the 
difference between both approaches is found to appear at the fifth order. Calculating the 
series up to the 200th order, the series still happens to be convergent, and no sign of a 
divergent behaviour is observed. As has already been discussed in [26], for k < 1/V, which 
corresponds to a s (M 2 ) < 2/j3\ = 0.44, the expansion ([4.1|) converges on the whole unit 



circle. Thus, also the FO series for R^'J® should be convergent for the physical value of 
a s (M 2 ), at least if only a finite number of c n \ coefficients are included. 

5. The general case 

In order to tackle the general problem with all /3-function coefficients f3\ to f3 4 being unequal 
to zero, the general expansion of a(— sqx) in terms of a So is employed, which takes the form: 



a(-s x) 



i + EE^M-*) 

i=l j=l 



(5.1) 



The coefficients dij can be calculated from the RGE for the QCD coupling eq. (UA). The 
coefficients of the highest and the lowest power of the logarithm can easily be given for ar- 
bitrary order and read d nn = (— (3i/2) n and d n \ = —f3 n /2. The remaining three coefficients 
up to fourth order are found to be: 

d 32 = | Pi , d 43 = - if 0{ 2 , d 42 = § (Pi + 20 x /? 3 ) . (5.2) 

In principle - given the necessary computing resources - the coefficients d%j can be cal- 
culated to an arbitrary order, taking into account the known /3-function coefficients, and 
setting the even higher-order ones to zero. 

The numerical analysis of the case including all available /3-function coefficients pro- 
ceeds along similar lines as the simplified example in the previous section. For CIPT, 
now the integrals of eq. ( |3.10j ) have to be calculated numerically by plugging a numerical 
solution of the RGE ( [2.4 ) for the running coupling and then performing the contour inte- 
gration. The FO result can be taken from eq. Q3.7| ). In addition including the fifth-order 
term, one obtains: 



D,CI 



R T 



r>D,FO 



a a 1 a 2 a 3 a A a 5 

1 + 0.1479 + 0.0297 + 0.0122 + 0.0047 + 0.0019 



3.5893, (5.3) 



1 + 0.1082 + 0.0609 + 0.0334 + 0.0144 + 0.0021 = 3.6571 . (5.4) 
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9 


10 


11 


12 


13 


14 


15 


16 


-0.004380 


0.005429 


0.010012 


0.006522 


-0.001223 


-0.006793 


-0.006409 


-0.001148 


17 


18 


19 


20 


21 


22 


23 


24 


0.004472 


0.005997 


0.002420 


-0.003022 


-0.005577 


-0.003107 


0.002045 


0.005211 


25 


26 


27 


28 


29 


30 


31 


32 


0.003582 


-0.001304 


-0.004970 


-0.004026 


0.000713 


0.004875 


0.004498 


-0.000213 



Table 2: Perturbative order n and corresponding terms in the expansion of R®' F /3. 



Again, the CI series displays a better convergence although the last included term in both 
expansions is almost of the same size. This time, however, the difference between CI and 
FO perturbation theory looks much more dramatic, as it amounts to 3 • 0.0226, more than 
ten times the last included term! This apparent ambiguity in the perturbative prediction of 
the hadronic t decay rate at the moment represents the dominant theoretical uncertainty 
in the extraction of the strong coupling a s from this channel. 

To gain further insight into the origin of the problem, like in the previous, simplified 
example, higher order contributions should be inspected. These can be acquired by ex- 
ploiting the expansion ( |5.1| ) for the running coupling in the contour integrals of eq. ( 3.1Q ) 
which can be expressed in terms of the Ji, and then reexpanding the result in terms of 
a S() . 3 The result of this exercise is presented in table ^. Like in the last section, after the 
first five positive terms, there are four negative ones after which the terms change sign at 



5 Of course, this is just another way of calculating higher order relations for the c n k coefficients, similar 



to the ones given in eq. (2.6) 
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Figure 2: Difference (R^ ,FO — R®' )/3 as a function of the perturbative order n up to which the 
terms in the FOPT series are summed. 
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every third order, and, all in all, decrease in magnitude. The negative sign of the sixth 
order term, which to some extent already lowers the difference between CIPT and FOPT, 
has also been observed recently in ref. [5]. 

Restricting oneself to the five c n \ coefficients discussed above as well as the four known 
/3-function coefficients, again the CI value ( |5,3D represents the full answer. This time, in 
figure H the difference between the FO expansion summed up to a particular order n and 
the complete CI result is displayed up to the 51st order. Until about the 30th order, the 
series shows the same behaviour as in the simplified example above. It oscillates around 
the CI result, while approaching it. However, for the general case the convergence to the 
CI result is much slower, and roughly after the 30th order, an asymptotic behaviour of 
the series sets in and the terms start to again increase. The observed asymptotic nature 
of the FOPT series for R® ,FO has nothing to do with the expected asymptotic behaviour 
of the c n i coefficients which is due to renormalons [37,38]. Inspection of the expansion 
( |5.1| ) reveals that for the used physical value of a s (M^), it diverges on part of the unit 
circle, leading to the behaviour detected in figure ||. The minimal deviation of the FO 
result at the 30th order is about 0.005, and the oscillation frequency is approximately six 
perturbative orders. Alas, again the maximum of the deviation occurs at the fifth order. 



6. Scale variations 



In perturbative QCD calculations, it is quite standard to estimate uncertainties due to as 
yet uncomputed higher orders by a variation of the in principle arbitrary renormalisation 
scale. Thus, in this section the influence of scale variations on FO and CI perturbation 
theory for the hadronic r decay rate shall be investigated. Let us begin with a study of 
renormalisation group improved FOPT. 

Like in ref. [26], the arbitrary renormalisation scale will be introduced by substituting 
// 2 = £ 2 so in the last line of eq. (|3.2| ), which yields: 

OO 71+1 k— 1 /, -. \ 

R?' FO (eso) = tfcE°"(* 2 *°)E***E( i jJii-^o^ 1 - 1 - (6.i) 

n=0 k=l 1=0 ^ ' 



The previous eq. (3J3) just corresponds to the special case £ = 1. One can then proceed 
with a numerical analysis, completely analogous to the presentation in the last section. 
The result of this exercise is shown in figure |3| up to the 60th perturbative order and for 
two values of £. The full triangles correspond to the case £ = 0.9 whereas the full circles 
result by setting £ = 1.1. To guide the eye, the data points have been connected by straight 
line segments. The required inputs for a(£ 2 M 2 ) have been calculated by solving the RGE 
(U) with the initial value a s (M 2 ) = 0.34. 

As can be observed from figure |3|, changing the renormalisation scale even by only a 
small amount, the behaviour of the FO series changes quite drastically. In the case of a 
smaller renormalisation scale £ = 0.9, the asymptotic behaviour already sets in around the 
15th order, and the minimal deviation from the CI result is as large as 0.014, while for the 
case of a larger renormalisation scale £ = 1.1, close inspection shows that the amplitude of 
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Figure 3: Difference (R^' FO — i?^' CI )/3 as a function of the perturbative order n up to which the 
terms in the FOPT series are summed and for an arbitrary renormalisation scale parameter £. Full 
triangles correspond to £ — 0.9; full circles to £ = 1.1. 



the oscillations decreases until about the 45th order after which it gradually again increases. 
This time the minimal amplitude turns out to be 0.0013. Like in the case £ = 1, however, 
for both examples the oscillation frequency is about six perturbative orders. Therefore, 
again for values of the renormalisation scale different from £ = 1, the FOPT series is found 
to be compatible with the CIPT value (|5.3|) , although due to its rather bad behaviour it 
should not be utilised for any phenomenological analysis. 

For CIPT, introducing an arbitrary renormalisation scale parameter £ can be achieved 
with the choice /j? = —^sqx in the first line of (p.2|), which results in 



n+l 



R?' C1 (es ) = AT c ^^(£ 2 S0 )E fcc »fe(-21n£) 



fc-i 



(6.2) 



n=0 



k=l 



This time, the equation has been rewritten in terms of the contour integrals </n(£ 2 so)> given 
by: 

1 f dr 

JZifso) = — j — (1 - xf (1 + x) a n (-fs x) . (6.3) 
l*l=i 

From eq. Q6.2j ) it is apparent that in contrast to CIPT with £ = 1 of eq. ( |3.9D , and similarly 
to FOPT, at a given order n, via the coefficients c n k, all independent coefficients en with i < 
n contribute. Even if again only the first five c n \ are assumed to be non-zero, nevertheless, 
the CIPT series for general £ does not terminate at this order. Hence, now also for CIPT 
one can investigate the difference between i?^ > ' CI (£ 2 so) for a given £ and the complete result 
R?' C1 of eq. Q. 

The result of this exercise is displayed in figure ||| up to the 20th perturbative order. 
The full triangles now correspond to the case £ = 0.5 whereas the full circles result by 
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Perturbative order n 

Figure 4: Difference (R^ ,ci (^ 2 sq) — R® ,cl )/3 as a function of the perturbative order n up to which 
the terms in the CIPT series are summed and for an arbitrary renormalisation scale parameter £. 
Full triangles correspond to £ = 0.5; full circles to £ = 2.0. 



setting £ = 2.0. To guide the eye, the data points have been connected by straight line 
segments. Generally, the CIPT for arbitrary £ is much more stable against a variation of 
the scale parameter than the corresponding FO series. Thus, a larger variation of £ could be 
considered. In the case of the smaller £ = 0.5, surprisingly, and most probably by chance, 
between the second and the fifth order, eq. fl6.2| ) already gives a rather good approximation 
to the full answer. Beginning with the sixth order, the deviation again fluctuates more 



strongly until it finally converges towards the result of eq. (5.3) for even higher orders. In 
the case of a larger scale with £ = 2.0, the approach towards the complete result is much 
more smooth. Coming from below, the difference only assumes small positive values before 
again approaching zero. Nevertheless, like for utilising FOPT, there is no real incentive 
to employ CIPT with values for £ different from one. This procedure only reshuffles the 
perturbative expansion, thereby transferring known contributions to higher orders which 
can as well be resummed. 

After the discussion of renormalisation scale dependence of FOPT and CIPT, let us now 
address the important question which are the uncertainties of the result ( |5.3|) and how they 
might be estimated. Obviously, one source of uncertainty resides in the contribution of yet 
uncalculated higher order coefficients c n \. Schemes to estimate these coefficients have been 
discussed in the literature [33-36] (and references therein) , 4 but it seems fair to say that no 
rigorous assessment of the corresponding uncertainty exists, prior to an explicit calculation 
of a certain coefficient. It also remains unclear if using the size of the last computed term 
as an estimation for yet higher order contributions is justified - although this is often done 
in practice. On the other hand, a variation of the renormalisation scale only amounts to 
a reordering of the perturbative expansion for a given set of known c n \ coefficients and 



4 Such estimates for the coefficients C41 and C51 had been employed above. 
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has nothing to do with the uncertainty due to so far uncalculated coefficients. This is 
somewhat different for a change of the renormalisation scheme, as will be discussed in the 
next section. 

A second source of uncertainty in CIPT results from the fact that the running of the 
QCD coupling is required while integrating over the contour in the complex s-plane, and 
the perturbative expansion for the /?-function has been truncated at the fourth order. An 
indication of the corresponding uncertainty can be gained by comparing the result ( |5.3| ) 
with the value that would be obtained while using three-loop running for a s . Numerically, 
the difference between both approaches is found to be 3 • 0.00099, only about half of the 
last included term in the perturbative series. The third, and final uncertainty on the result 
( |5.3| ) comes from the error of the input value for a s (M^). This uncertainty, however, can 
be worked out straightforwardly, assuming a given error on a s (M^). 



7. Scheme variations 

Besides investigating the dependence of the perturbative expansion on the choice of the 
renormalisation scale, also a modification of the renormalisation scheme provides interesting 
further information. In analogy to the general expression Q5.1| ) for the shift of the scale at 
which the QCD coupling is evaluated, a transformation of the renormalisation scheme can 
be represented by the following equation: 



1 + di 



a s 



(7.1) 



As before, initial quantities, like the coupling a So , are given in the MS-scheme, whereas 
tilded quantities will always correspond to the transformed scheme. For the case of the 
modification of the renormalisation scheme, the same scale so is assumed in both couplings 
a So and a SQ . In principle, one could have considered both transformations simultaneously, 
in which case the coefficients di would have corresponded to coefficients dio in eq. ([xl]). 
but it appears more transparent to discuss them separately. 

Because the Adler function is a physical quantity, it should not depend on the renor- 
malisation scheme used. This can only be true, if together with the coupling also the 
coefficients of the perturbative expansion c n \ get modified. Comparing in eq. fl2.2| ) or 
(3.10) the expansion in terms of the two schemes, it is straightforward to read off the 



transformed coefficients [34]. The scheme dependence of the coefficients only starts at or- 
der a 2 , which implies coi = coi and c\\ = c\\- Up to order a 5 the remaining relations are 
found to be: 

C21 = c 2 i - d\c\\ , c 3 i = c 3 i - 2dic 2 i + (2d\ - d 2 ) c n , 
c 4 i = c 4 i - 3ffic 3 i + (5df - 2d 2 ) c 2 i - (5cZf - 5d\d 2 + d 3 ) c\\ , (7.2) 
^5i = C5i - 4diC4i + {9dj - 3d 2 ) c 3 i - (14d? - 12cM 2 + 2d 3 ) c 2 i 
+ (14^ - 21dld 2 + 3dj + 6cW 3 - d A ) c u . 
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This shows that, even though CIPT still takes the functional form of eq. (|3,9| ), now ex- 
pressed in terms of the new coefficients c n \ as well as new contour integrals J%(so), assum- 
ing the series to terminate in one scheme, this is no longer the case in a different scheme. 
Through relations of the type ( |7.2[ ), coefficients c n \ are generated to all orders. However, 
to be able to tell in which scheme a certain coefficient c n \ vanishes, it necessarily has to be 
calculated in at least one reference scheme. 

Like for eq. ( |7.2[ ) , it is a simple matter to derive from ( |7.1[ ) the RGE for the coupling a So 
and the corresponding (3- function coefficients fy. Thereby, one can confirm the well-known 
fact that the first two coefficients of the /^-function are independent of the renormalisation 
scheme, i.e. (3\ = j3\ and $ 2 = /?2- 5 Starting with j3\, it is found that the /3 coefficients 
depend on the di, such that the coefficient (3^ depends on all di with i = f, . . . ,k — 1. 
Explicitly, the relations for (3% and (3$ take the form: 

fa = A - dift - (4 - <h)Pi , fa = At - 2di/? 3 + d\(3 2 + (44 - 6cW 2 + 2d 3 )A • (7.3) 

Again, one could define a scheme, in which all (3k = for k > 3 [40]. But in order to do 
this, one needs to know the /3-function coefficients in one reference scheme. 

Let us now proceed with a numerical analysis of the scheme dependence. The principle 
of minimal sensitivity (PMS) [33] attaches a special importance to the point where the 
derivatives of the physical quantity with respect to the scheme-parameters di vanish, that 
is where the physical quantity, with respect to the di has an extremum or a saddle point and 

D CI 

thus the scheme dependence is weakest. Performing this exercise for R T ' as a function 
of d\ and d 2 , and putting the remaining di to zero, one finds a maximum at d\ = 0.165 
and d 2 = 0.210. (Trying to also include d^, in the analysis, it is found that the extremum 
drifts off to values of the d-coefficients, where the scheme transformation ([O]) is no longer 
perturbative.) 

In figure |j| the dependence of R T ' , including up to the fifth perturbative order, on 
the scheme parameters d\, d 2 and d% is displayed in a graphical form. To this end, always 
only one parameter is varied and the remaining parameters are kept fixed either at the 
maximum given above for d\ and d 2 or at zero for d% and d±. In order that the scheme 
transformation stays perturbative, a reasonable choice for the variation of the di seems to 
be the range — 1 < di < 1, such that the correction is always at most about 10% of the 
previous term. The variation of d\ is given by the solid line, with the maximum being at 
d\ = 0.165. Next, the variation of d 2 is given by the dashed line, with the maximum at 
d 2 = 0.210, and finally the variation of cfe corresponds to the dotted line. 

Several remarks are in order. As expected, for higher corrections, the dependence on 
the scheme parameters di gets weaker and weaker, being strongest for d\ and weakest for 

D CI 

c/3 . Actually, varying d%, the change in R T ' is only by one unit in the fifth digit. While 



varying d\, the maximal difference to the result of eq. (5.3) is given by 3 • 0.00073. Thus, 
the scheme variation is of a similar size as the expected uncertainty from higher orders in 
the /3-function, considered in the last section. Furthermore, the "optimal" scheme, where 
the dependence on the parameters di is weakest lies rather close to the MS scheme in which 



At least in renormalisation schemes which are mass and gauge-parameter independent [39]. 
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Figure 5: R®' Cl as a function of the renormalisation scheme parameters d\, c?2 and d%. The solid 
line corresponds to a variation of d± , the dashed line to di and the dotted line to dj, . 



all di = 0. Therefore, also the speed of convergence of the CI perturbative series in the MS 
scheme can be expected to be close to optimal. 



8. Conclusions 

Until today, the extraction of a s from the hadronic r decay rate is hampered by an apparent 
ambiguity between performing the renormalisation group improvement of the perturbative 
series at a fixed order or in the so-called contour-improved scheme [25,26]. This ambiguity 
represents the dominant theoretical uncertainty for a s as extracted from hadronic r decays. 

Further insight into the origin of the problem can be gained by pursuing the fixed-order 
expansion to larger orders, still staying consistent with the terms that have been resummed 
in contour-improved perturbation theory. In a simplified example, where only the first 
coefficient of the /3-function, fix, is kept non- vanishing, this could be done analytically 
in closed form. For the more general case including the four currently known terms of 
the /3-function, the expansion of the fixed-order result can be calculated to, in principle, 
arbitrarily high order. The behaviour of the fixed order series in the simple example as well 
as in the more general case discussed in sections 4 and 5 have been presented in tables [j] 
and H respectively, and graphically illustrated in figures || as well as |2[ 

It is observed that in the full QCD case the fixed-order result oscillates around the 
contour-improved value with an oscillation frequency of approximately six perturbative or- 
ders, converging to it until the 30th order after which an asymptotic behaviour of the series 
shows up and the terms again increase. Just including five perturbative terms in the fixed- 
order case happens to turn out most unfortunate, as for this order the deviation between 
contour-improved and fixed-order results has a maximum. Analysing the comparison of 
fixed-order as well as contour-improved perturbation theory for an arbitrary renormalisa- 
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tion scale parameter £ unequal to one reveals that also in this case both are compatible. 
However, invoking £ 7^ 1 in CIPT reshuffles known contributions to higher orders thus 
throwing away available information. Changing the renormalisation scheme is a somewhat 
different issue, for there the independent coefficients c n \ themselves get modified. Still, 
this does not allow for an unambiguous assessment on the uncertainties resulting from so 
far neglected higher orders. 

All the discussion above demonstrates, that the renormalisation group improvement of 
perturbative series in finite-energy sum rules, like for the hadronic r decay rate in particular, 
should be performed in the contour-improved scheme. This approach also is most natural 
in the sense that all terms proportional to a particular perturbative coefficient c n \ in the 
correlation function are resummed to yield the contour-improved contribution at order n, 
whereas in the fixed-order approach terms proportional to the coefficient c n \ appear at 
all orders equal or greater to n. A discussion of renormalisation group improvement of 
quark-mass corrections in finite-energy sum rules along similar lines will be presented in a 
forthcoming publication. 
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